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This school is not about N-body simulations. 
But rather, about the art of N-body simulations. 


To run an N-body simulation you require a working code and a computer. 


To do an N-body simulation and actually learn something from it, you 
also need understanding of the physics and the numerics involved. 


After all, computers only give you numbers, but not insight. 
Running simulations can actually lead you on the wrong path. 


It is only by doing simulations that you can gain understanding. 


It is the goal of this school to show you how to do N-body simulations. 
This is an art that, unfortunately, is usually not taught in books or 
courses, but it is passed from tutors to students. 
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Background 


To talk about N-body simulations, we should go back to gravity. 


It all began with an apple (or so thay say) and Newton in the XVII century. 


“Planetas omnes in se muto graves esse jam ante probavimus, 
ut & gravitatem in unumquemque seorsim spectatum 

esse reciproce ut quadratum distantiœ locorum a centro planetæ. 
Et inde consequens est gravitatem in omnes proportionalem 
esse materiæ in iisdem”. 


Isaac Newton, 
Propositio vii, Teorema VII. 
Philosophiæ Naturalis Principia Mathematica. 
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Background 


To talk about N-body simulations, we should go back to gravity. 


The solution of the two-body problem, lead Newton to the invention of Calculus. 
The three-body problem, however, proved intractable. Newton wrote that its solution, 


“Unless | am much mistaken, it would exceed the force of human 
wit to consider so many causes of motion at the same time’. 


And right there, from the very beginning, laid the seed that would lead to N- 
body simulations. 
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Background 


The historical development of the field of classical gravitation 


The gravitational N-body problem (in the classical regime), has a historical 
development that is inverted with respect to other branches of Physics: the 
rest has started from a global description of phenomena to develop, 
subsequently, a local description of the fundamental interaction among the 
components of the system. 


Thermodynamics œ Statistical Mechanics, 
Fluid Mechanics = Kinetic Theory of gases, 
Bulk properties of materials => Atomic and molecular structure. 
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Background 


The historical development of the field of classical gravitation 


In the gravitational problem we knew the interaction rule | 
from the beginning. TENGE 


The difficulty was in applying it to several particles at 
once. 


The study of the N-body problem followed two paths: the first relied on 
finding approximations to the true solution. This gave rise to Celestial 
Mechanics. The second approach looked for solutions transforming the 
original problem in another, equivalent one. This developed into Analytic 
Mechanics or Rational Mechanics. 

In the last 100 years, the study of the N-body problem lead to the 
development of another area of inquiry: Stellar Dynamics. 
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Background 


e Pierre Simon de Laplace investigated the stability of the Solar System 
and published in 1799 his classic work: Mécanique Celeste. 


e Leonhard Euler seems to have been the first to study the three-body 
problem as a pure mathematical problem. He found a solution of the 
restricted 3-body problem in which all particles move along a straight 
path. 


e Joseph Louis Lagrange studied the motion of the Moon and gave a 
detailed analysis of the orbits of Jupiter's satellites. He confirmed Euler's 
solution for the general 3-body problem and found an additional solution 
where the 3 bodies move forming an equilateral triangle. 


e Karl Friedrich Gauss developed a general method to determine orbits 
based on only three observations and found the asteroid Ceres, which 
had been lost when it passed behind the Sun. 


e Henri Jules Poincare studied very thoroughly the 3-body problem. The 
methods he introduced to study its stability lead to the development of 
Topology and to the discovery of the concept of Chaos. 
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Background 


The N-body problem is of great interest for several sciences 


(i) M = T*W with canonical symplectic structure, 


W=R°xR*\A, A= {(¢q) lq € R°} 


(ii) m € M (initial conditions); 
(iii) u € R, u > 0 (mass ratio); 


(iv) H” € F(M) given by 
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Background 


The N-body problem and the quest for power among European powers 


The great European powers established national 
observatories and promoted the development of 
mathematics and astronomy, trying to solve a 
problem of large strategic importance: navigation. 


The search for a solution split in two camps: 

Some turned to mechanical devices to measure 
time accurately and thus determine geographical 
position, while others kept looking at the sky for the 
perfect clock. 
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Background 


The N-body problem and the quest for power among European powers 


The second approach motivated the study of the N-body problem and brought 
an enormous development of Mathematics and the Physical Sciences. 


International School on Galactic and Cosmological N-Body Simulations i 
July 23 — August 5, 2006 Tonantzintla — Puebla, Mexico The Art oF N-Body 


Luis Aguilar 


Background 


So, it is fair to say that the N-body problem is one of the 
most fruitful problems in the history of science. 
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The nature of gravity 


Although gravity does not have an intrinsic length scale, its application 
to discrete systems does introduce a characteristic scale. 
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The nature of gravity 


Although gravity does not have an intrinsic length scale, its application 
to discrete systems does introduce a characteristic scale. 


Let's imagine a test particle at the center of a spherical system that contains a huge 
number of point particles. 


It is clear that the net average force on the test particle will be zero, since the 
distribution of matter is symmetric around it: there is, on average, the same number of 
particles pulling on each direction. 
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The nature of gravity 


Although gravity does not have an intrinsic length scale, its application 
to discrete systems does introduce a characteristic scale. 


Let's imagine a test particle at the center of a spherical system that contains a huge 
number of point particles. 


It is clear that the net average force on the test particle will be zero, since the 
distribution of matter is symmetric around it: there is, on average, the same number of 
particles pulling on each direction. 


The instantaneous force, however, will fluctuate 
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The nature of gravity 


Although gravity does not have an intrinsic length scale, its application 
to discrete systems does introduce a characteristic scale. 


If we center our attention to a particular region which contains N stars at a distance r 
from us, it is clear that the force due to it is proportional to: Fœ N 
2 


Due to the finite discrete nature of the system, 
there will be fluctuations in the force of a OF œ oN = 
magnitude given by: g 


JN 


r2 
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The nature of gravity 


Although gravity does not have an intrinsic length scale, its application 
to discrete systems does introduce a characteristic scale. 


If we center our attention to a particular region which contains N stars at a distance r 
from us, it is clear that the force due to it is proportional to: F œ N 
2 
3 


Due to the finite discrete nature of the system, T 
there will be fluctuations in the force of a OF «2= — 
magnitude given by: a r 


If we now look at another region at a distance 1’ 
and ask that it exerts the same average force, it 
is clear that it must contain a number of stars 

given by: 


Nice r 
and its fluctuations will be given by, 
ye eee 
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AF Be 1 This is a very important result, it means that if we split a discrete 
r system into parcels, such that each one exerts on average the same 
force on the test particle, the fluctuations in their forces will shrink wit! 
distance as 1/r. 
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The nature of gravity 


AF Be 1 This is avery important result, it means that if we split a discrete 
r system into parcels, such that each one exerts on average the same 
force on the test particle, the fluctuations in their forces will shrink wit! 
distance as 1/r. 


We can then split the gravitational influence into two asymptotic 
regimes: a fluctuating part due to close neighbors and a smooth part 
due to distant particles. 


irregular smooth 
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The nature of gravity 


The characteristic time in which the irregular part dominates is called the 
collisional time and it can be shown that it is of the form: 


ta © (Æ) N /log(N) 


ay 


On the other hand, the natural time unit in a gravitational system is the so 
called dynamical time in which particles cross the system: 


R 
cc ——_— 
l ayn ( <v? > ) 


Their ratio measures the importance of the irregular part with respect to 
the smooth one: 


t/t, œN /log(N) 


col dyn 
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The nature of gravity 


The smaller this ratio, the more the system will be dominated by the irregular regim: 
where the behavior of the force is dominated by collisions with close neighbors. Tht 
we can split self-gravitating N-body systems in two types: collisional and collisionle: 
systems; the fundamental difference between them being set by the number of 


particles. 


to (10° years) ln (10° years) 
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The nature of gravity 


The fact that the importance of collisions decreases as N increases is a bit 
counterintuitive, so it deserves a bit of thought 


Let us imagine a self-gravitating system of N particles in virial equilibrium, it is clear 
that as time passes some collisions will occur. 
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The nature of gravity 


The fact that the importance of collisions decreases as N increases is a bit 
counterintuitive, so it deserves a bit of thought 


Let us now imagine that we split each particle in two, thus doubling the 
number of particles in the system without perturbing the virial equilibrium. 
It is clear that the number of collisions per unit time will double. 


So, how can we say 
that the importance of 
collisions diminishes? 
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The nature of gravity 


What matters is not the number of collisions, but the effect they have. 


In the case of cars, collisions imply physical contact, and the faster the relative 
velocity, the larger the damage produced. 


> + 4 


In the case of stars, we have long range interactions, where the larger the 
relative velocity, the shorter the interaction time, and less damage will occur. 


> + 


fr fr 
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The nature of gravity 


In the case of our hypothetical system with N particles, the strength of the gravitational 
force goes as the product of the masses divided by their separation: 


F «<(mem)/r? 


When we split each particle in two, we doubled the number of particles (and collisions), bi 
the mass was halved, so the numerator got reduced by a factor of 1/4: 


m>m/2—> FO F/2° 


Its true that the mean separation between particles got reduced as well, which drives the 
effect in the other direction, but this effect can’t overcome the mass effect: 


re NU? = F+>Fx(1/27)(2?°?)=Fx27°° 


Even after multiplying by 2 to account for the larger number of collisions, the net effect shrir 
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The basic parts of an N-body code 


The two most important parts of an N-Body code are: 


e A force routine that computes forces. 


e An integrator that integrates the equations of motion. 


The choices made for these two parts largely define the merits and 
limits of an N-body code. 
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The basic parts of an N-body code 


Force calculators 


The most important feature that defines the type of force calculator is 
the collisionality of the system. Other properties, like overall symmetry, 
are important too. 


Less collisional —— — — — 


PP PSM Tree MFC PM Fourier 
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The basic parts of an N-body code 


Integrators 


Integrators are defined to a large extent by the numerical precision 
required. The dynamical range of the timescales as well as restrictions 
on computer memory or computing speed are also relevant. 


Fixed time-step vs variable time step 
| Low order vs high order schemes 
_ Symplectic vs conventional integrators 
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The limits of an N-body simulation 


We must never forget that a given N-body simulation, however large and 
carefully planned, is always an approximation: 


| e Limitations of the integrator 
e Approximations of the force calculator 
e Finite dynamical range 
e Issues related to coarseness of the simulation 
compared to real systems 
e Non-equilibrium initial conditions 
_ e Convergence 


You would be surprised to see the number of published works that are 
flawed, not for limitations of computer technology, but by poor design. 


International School on Galactic and Cosmological N-Body Simulations 
July 23 — August 5, 2006 Tonanizintla — Puebla, Mexico 


The Art of N-Body 
z eel 


Instituto de Astronomia UNAM E 


Outline of the school 


Topics and lecturers: 


Luis Aguilar 


e Unix & Nemo (Jean Charles Lambert & Peter Teuben) 
e Initial conditions (Luis Aguilar) 

e Force routines (Walter Dehnen) 

«Integrators (Peter Teuben) 

e The Astrophysics behind (Lia Athanassoula) 

e Cosmology (Julio Navarro) 


Walter Dehnen Julio Navarro 


Lia Athanassoula 
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Outline of the school 


Topics and lecturers: 


e Unix & Nemo (Jean Charles Lambert & Peter Teuben) 
e Initial conditions (Luis Aguilar) 

e Force routines (Walter Dehnen) 

«Integrators (Peter Teuben) 

e The Astrophysics behind (Lia Athanassoula) 

e Cosmology (Julio Navarro) 


Structure of the school: 


e Regular morning lectures 
e Afternoon computer labs 
e Evening conferences: 
- A. Bosma, |. Puerari, O. Valenzuela, P. Colin, V. Avila 


Hands on experience: 


e Basic setup 
e Simple exercises 
eSchool projects 
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My favorite N-body simulation 


| can not finish without telling you about my favorite simulation 


The year is 1941 and based on surveys like the Shapley-Ames, people were 
aware of the tendency of luminous galaxies to gather in large conglomerates 
called clusters of galaxies. 


Erik Holmberg conjectured that energy dissipation during galaxy-galaxy 
encounters was responsible, at least in part, for this tendency. 


However, he needed to perform a simulation to gauge the magnitude of the 
effect. Since at this time no computers were available, at least not to him, he 
had to do it some other way. 
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My favorite N-body simulation 


Holmberg replaced gravity by light to perform the first published “N-body 
simulation” in Galactic Dynamics. 


THE ASTROPHYSICAL JOURNAL 


AN INTERNATIONAL REVIEW OF SPECTROSCOPY AND 
ASTRONOMICAL PHYSICS 


VOLUME 94 , NOVEMBER 1941 NUMBER 3 


ON THE CLUSTERING TENDENCIES AMONG THE NEBULAE 


II. A STUDY OF ENCOUNTERS BETWEEN LABORATORY MODELS OF 
STELLAR SYSTEMS BY A NEW INTEGRATION PROCEDURE 


ERIK HOLMBERG 


ABSTRACT 


In a previous paper! the writer discussed the possibility of explaining the observed clustering effects 
among extragalactic nebulae as a result of captures. The present investigation deals with the important 
problem of whether the loss of energy resulting from the tidal disturbances at a close encounter between 
two nebulae is large enough to effect a capture. The tidal deformations of two models of stellar systems, 

ssing each other at a small distance, are studied by reconstructing, piece by piece, the orbits described 
y the individual mass elements, The difficulty of integrating the total gravitational force acting upon 
a certain element at a certain point of time is solved by replacing gravitation by light. The mass elements 
are represented by light-bulbs, the candle power being proportional to mass, and the total light is meas- 
ured by a photocell (Fig. 1). The nebulae are assumed to have a flattened shape, and each is represented 
by 37 light-bulbs. It is found that the tidal deformations cause an increase in the attraction between the 
two objects, the increase reaching its maximum value when the nebulae are separating, i.e., after the 
passage. The resulting loss of energy (Fig. 6) is comparatively large and may, in favorable cases, effect 
a capture. The spiral arms developing during the encounter (Figs. 4) represent an interesting by-product 
of the investigation. The direction of the arms depends on the direction of rotation of the nebulae with 
respect to the direction of their space motions. 
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The Ead 


or rather, the beginning. 
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